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Abstract 

We apply the equal load-sharing fiber bundle model of fracture failure in composite 
materials to model the traffic failure in a system of parallel road network in a city. 
For some special distributions of traffic handling capacities (thresholds) of the roads, 
the critical behavior of the jamming transition can be studied analytically. This has 
been compared with that for the asymmetric simple exclusion process in a single 
channel or road. 
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1 Introduction 



Traffic jams or congestions occur essentially due to tlie excluded volume effects 
(among the vehicles) in a single road and due to the cooperative (traffic) load 
sharing by the (free) lanes or roads in multiply connected road networks (see 
e.g., [1]). We will discuss here briefly how the Fiber Bundle Model (FBM) 
for fracture failure in composite materials (see e.g, [2]) can be easily adopted 
for the study of (global) traffic jam in a city network of roads. In the equal- 
load-sharing FBM, the load carried by the failed fibers is transferred and 
shared uniformly by the other surviving fibers in the bundle and therefore 
the dynamics of failure propagation in the system is analytically tractable [2]. 
Using this model for the traffic network, it is shown here that the generic 
equation for the approach of the jamming transition in FBM corresponds to 
that for the Asymmetric Simple Exclusion Processes (ASEP) leading to the 
transport failure transition in a single channel or lane (see e.g., [3]). 
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Fig. 1. The highway traffic current Iq gets fragmented into i's and the narrower 
roads having threshold current Zc„ < i gets congested or blocked and this extra 
traffic load per uncongested roads get equally distributed, causing further blocking 
of some more roads. 



2 Model 



Let the suburban highway traffic, while entering the city, get fragmented 
equally through the various narrower streets within the city and get combined 

again outside the city (see Fig. 1). Let the input traffic current be denoted by 
Jo and the total output traffic current be denoted by It- In the steady state, 
without any global traffic jam. It = lo- In case. It falls below /q, the global 
jam starts and soon It drops to zero. This occurs if Iq > the traffic cur- 
rent in the network beyond which global traffic jam occurs. The steady state 
flow picture then looks as shown in Fig. 2. Let us assume that the parallel 
roads within the city have got different thresholds for traffic handling capac- 
ity: ic^,ic2, ■ ■ ■ , icN the N different roads (the n-th road gets jammed if the 
traffic current i per road exceeds ic„). Initially i = Iq/N and increases as some 
of the roads get jammed and the same traffic load 7, has to be shared equally 
by a lower number of unjammed roads. Next, we assume that the distribution 
p{ic) of these thresholds is uniform (see Fig. 3) upto a maximum threshold 
current (corresponding to the widest road traffic current capacity), which is 
normalized to unity (sets the scale for 1^) . 

The jamming dynamics in this model starts from the n-th road (say in the 
morning) when the traffic load i per city roads exceeds the threshold ic,, of that 
road. Due to this jam, total number of intact or uncongested roads decreases 
and rest of these roads have to bear the entire traffic load in the system. 
Hence effective traffic load or stress on the uncongested roads increases and 
this compels some more roads to fail or get jammed. These two sequential 
operations, namely the stress or traffic load redistribution and further failure 
in service of roads continue till an equilibrium is reached, where either the 
surviving roads are strong (wide) enough to share equally and carry the entire 
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Fig. 2. Traffic jam (global) in the network of roads occurs when Iq > Ic, the 
threshold traffic capacity of the city road-network. 

traffic load on the system (for Iq < Ic) or all the roads fail (for Iq > Ic) and a 
(global) traffic jam occurs in the entire) road network system (to be brought 
back perhaps in the night when the load Iq decreases much below Ic)- 



3 Jamming dynamics 



This jamming dynamics can be represented by recursion relations in discrete 
time steps. Let us define Ut{i) to be the fraction of uncontested roads in the 
network that survive after (discrete) time step t, counted from the time t — 
when the load (at the level Iq = iN) is put in the system (time step indicates 
the number of stress redistributions). As such, Ut{i = 0) = 1 for all t and 
Ut{i) = 1 for t = for any i; Ut{i) = U*{i) ^ for t — oo if /o < h, and 
Ut{i) = for t ^ oo if /o > Ic- 

Therefore Ut{i) follows a simple recursion relation (see Fig. 3) 
TT ^ . . lo 

or, = A. (1) 

In equilibrium Ut+i = Ut = U* and thus (1) is quadratic in U* : 
U*' -U* + i^ 0. 



The solution is 

U*{i)^l±{ic-i)'/';ic^\. 
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Fig. 3. The simple model considered here assumes uniform density p{i) of the road 
traffic flow capacity distribution up to a cutoff capacity (normalized to unity). At 
any traffic load per road level it at time t, the fraction it of roads fails and fraction 
1 — it survives. 

Here ic — h/N is the critical value of traffic current (per road) beyond which 
the bundle fails completely. The quantity U*{i) must be real valued as it has a 
physical meaning: it is the fraction of the roads that remains in service under 
a fixed traffic load i when the traffic load per road lies in the range < i < ic- 
Clearly, ^7*(0) = 1. The solution with (+) sign is therefore the only physical 
one. Hence, the physical solution can be written as 



For i > we can not get a real- valued fixed point as the dynamics never stops 
until Ut — when the network gets completely jammed. 



4 Critical behavior 

(a) At i < ic 

It may be noted that the quantity U*{i) — U*{ic) behaves like an order pa- 
rameter that determines a transition from a state of partial failure of the 
system (i < to a state of total failure (i > ic) : 




(2) 



= U*{i)-U*{ic)^{ic-if](5 
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To study the dynamics away from criticality (i ic from below), we replace 
the recursion relation (1) by a differential equation 

_dU _ W^-U + i 
~~dt ~ U ■ 

Close to the fixed point we write Ut{i) — U*{i) + e (where e — >^ 0). This gives 

e = Ut{i)-U*{i)^eM-tlr), (4) 



where t = | \{ic — i) + 1 . Thus, near the critical point (for jamming 
transition) we can write 



T (X {ic — i) a = -. 



(5) 



Therefore the relaxation time diverges following a power-law as i ^ ic from 
below. 

One can also consider the breakdown susceptibility x, defined as the change 
of U*{i) due to an infinitesimal increment of the traffic stress i 



X 



dU*{i) 



di 



I,. 1 



(6) 



Hence the susceptibility diverges as the applied stress i approaches the critical 
value ic — \. Such a divergence in x can be clearly observed in the numerical 
study of these dynamics. 



(b) At I 



At the critical point {i 



we observe a different dynamic critical behavior 



in the relaxation of the failure process. From the recursion relation (1), it can 
be shown that decay of the fraction Ut{ic) of uncongested roads that remain 
in service at time t follows a simple power-law decay : 



[/, = 1(1 + ^), 
2^ t + V' 



(7) 



starting from Uq = 1. For large t {t ^ oo), this reduces to Ut — 1/2 (x t 
5 — 1; a. power law, and is a robust characterization of the critical state. 
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Fig. 4. The transport current I in the one dimensional lane or road is possible if, 
say, the i-th site is occupied and the {i + l)-th site is vacant. The inter-site hopping 
probability is indicated by J. 

5 Universality class of the model 



The universality class of the model can be checked (cf. [2]) taking two other 
types of road capacity distributions p{i): (1) linearly increasing density distri- 
bution and (11) linearly decreasing density distribution of the current thresh- 
olds within the limit and 1. One can show that while ic changes with different 
strength distributions {ic = y^4/27 for case (I) and ic = 4/27 for case 11), the 
critical behavior remains unchanged: a — 1/2 — (3 — j, 5 — 1 ior all these 
equal (traffic) load sharing models. 



6 Transport in ASEP 



Let us consider the simplest version of the asymmetric simple exclusion process 
transport in a chain (see Fig. 4). Here transport corresponds to movement of 
vehicles, which is possible only when a vehicle at site i, say, moves to the 
vacant site i + The transport current / is then given by (see e.g., [3]) 

7 = Ja(1-A+i), (8) 



where pi denotes the site occupation density at site i and J denotes the inter- 
site hopping probability. The above equation can be easily recast in the form 

= 1 - -, (9) 

Pi 



where a = I/J. Formally it is the same as the recursion relation (1) for the 
density of uncongested roads in the FBM model discussed above; the site index 
here in ASEP plays the role of time index in FBM. 
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7 Summciry and discussion 

We have applied here the equal load-sharing fiber bundle model of fracture 
failure in composite materials to model the traffic failure in a network of 
parallel roads in a city. For some special distributions of traffic handling ca- 
pacities (thresholds) of the roads, namely uniform and also linearly increasing 
or decreasing density of road capacities, the critical behavior of the jamming 
transition has been studied analytically. This dynamics of traffic jam progres- 
sion (given by the recursion relation (1), for uniform density of road capacity, 
as shown in Fig. 3 ) has been compared with that for the asymmetric simple 
exclusion process (9) in a single channel or road; the exact correspondence 
indicates identical critical behavior in both the cases. The same universality 
for other cases of p{i) in FBM suggests similar behavior for other equivalent 
ASEP cases as well. 
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